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I.  Introduction 

In  all  nonlinear  oscillators,  the  growth  of  the  amplitude  of  the  oscillation  is  ultimately 
curbed  by  back  reaction  of  the  growing  signal  on  the  agent  powering  the  oscillation.  In  the  case 
of  the  semiconductor  laser  the  back  reaction  consists  in  the  increased  recombination  rate  of  the 
excess  carriers  responsible  for  the  laser  action.  At  any  point  in  the  device,  the  increase  in 
recombination  rate  over  its  normal  equilibrium  value  is  proportional  to  the  local  power  density 
of  the  laser  field.  This  depletion  of  the  carrier  density  diminishes  the  gain  of  the  device, 
providing  the  negative  feedback  that  ultimately  stabilizes  the  laser  output  at  a  self-consistent 
operating  level.  Since  to  a  good  approximation  the  increase  in  recombination  rate  is  a  local 
effect,  the  excess  carrier  density  profile  will  be  nonuniform,  reflecting  the  intensity  distribution 
of  the  mode,  but  this  non-uniformity  is  counteracted  by  carrier  diffusion. 

Thus  one  has  a  nonlinear  problem  to  solve:  Maxwell’s  equations  with  carrier  density 
dependent  dielectric  constant,  together  with  the  diffusion  equation  with  field  intensity  dependent 
recombination  rate,  and  with  a  steady  source  term  reflecting  the  carrier  injection  process  in  the 
form  of  an  assigned  current.  One  can  go  back  further  and  trace  that  current  back  to  a  boundary 
value  problem  for  the  voltage  distribution  that  sets  up  that  current,*  making  extensive  use  of  the 

computer.  We  stop  short  of  this  step  since  it  is  difficult  to  see  how  it  can  affect  the  results  in  a 
major  way. 

The  method  we  use  to  "solve"  the  simultaneous  nonlinear  equations  is  conventional:  the 
electric  field  satisfying  Maxwell’s  equations  with  a  complex  dielectric  constant  (depending  on 
carrier  density)  and  the  carrier  density  are  expanded  in  a  series  of  cavity  modes,  with  amplitudes 
that  are  functions  of  time.  Then  the  condition  for  a  steady  state  is  written  down  (in  the  steady 

*H.  Suhl  is  Professor  of  Physics  at  the  University  of  California,  San  Diego.  He  was  supported  in  this  work  by 
USAF  Contract  F04701-88-C-0089. 


1 


state,  the  electric  field  varies  with  time  as  a  superposition  of  sinusoids,  whereas,  to  a  very  good 
approximation,  the  carrier  density  is  time  independent).  The  steady  state  condition  yields  a 
sequence  of  recursion  relations  for  the  Fourier  coefficients,  which  is  broken  off  at  an  early  stage. 
This  is  justified  if  the  field-induced  recombination  rate  (the  ultimate  source  of  the  nonlinearity)  is 
small  compared  with  the  equilibrium  recombination  rate. 

There  is  a  large  body  of  literature  on  this  subject.3-6  Most  of  it  addresses  specific  issues 
of  diode  laser  operation.  In  this  paper,  we  analyze  only  the  most  primitive  “toy  model”  that  still 
exhibits  the  spatial  hole  burning  phenomenon  in  analytic  (as  distinct  from  purely  computational) 
form. 


11.  Analysis 

Let  n  denote  the  excess  carrier  density.  We  assume  that  the  diffusion  equation  for  n  can 
be 


—  =  -  +  D^n  -  -  -  2aF^n  (1) 

dt  6  X 

where  J  is  the  driving  current  (in  general  a  function  of  position),  5  is  the  thickness  of  the 
active  layer,  a  is  a  coupling  constant  to  the  electromagnetic  field  and  D  is  the  diffusion 
coefficient.  (The  sign  of  7  is  as  shown,  because  carriers  are  being  generated  by  it.)  ~P~  is 
the  time  averaged  square  of  the  total  electromagnetic  field  (time  averaged,  because  the  carrier 
concentration  is  unlikely  to  follow  variations  at  optical  frequencies).  If  n  were  the  total 
carrier  density,  Eq.  (1)  would  have  to  include  higher  powers  of  n  inasmuch  as  the 
recombination  process  is  bimolecular  and  possibly  also  has  an  Auger  component  that  is  cubic 
in  However,  if  n  is  regarded  as  the  excess  carrier  density  needed  to  push  the  system  over 
the  threshold  to  its  operating  point,  and  that  excess  is  not  too  great,  equation  (1)  is  adequate. 
To  begin  with,  consider  a  cavity  consisting  of  two  infinite  parallel  plates  a  distance  a  apart, 
with  the  E-field  polarized  parallel  to  the  plates.  For  this  polarization  direction,  and  in  the 
absence  of  injected  carriers,  there  can  be  resonances  at  a  series  of  frequencies:  the  lowest 
is  at  w  =  r/((t/t(p.),  and  the  higher  ones  are  integral  multiples  thereof.  Here  €<,  is  the 
real  part  of  the  dielectric  constant  of  the  unactivated  medium.  In  the  presence  of  the  non- 
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linearity  implied  by  the  last  tenn  in  Eq.  (1),  the  steady  state  cannot  involve  one  single  mode 
alone.  We  shall  restrict  the  argument  to  a  superposition  of  two  modes.  In  the  presence  of 
the  current  J,  the  dielectric  constcint  becomes  £1,2  =  £0  +  +  t(7o  “  71.2’^)]  =  to 

first  order  in  the  excess  carrier  concentration  n.  Here  70  denotes  the  intrinsic  loss  in  the 
cavity,  and  71,2  the  gains  in  the  two  modes,  per  unit  carrier  concentration.  For  simplicity 
we  assume  that  the  real  part  of  the  dielectric  constant  is  the  same  for  both  modes.  The 
equation  for  the  Enfield  for  each  mode  is 

d^E  ^ 

dz^  ~  dt^  ’ 


where  2  is  the  distance  across  the  cavity.  For  real  e  it  has  a  solution  of  the  form  F(2)e’“*. 
When  e  is  complex,  but  its  imaginary  part  is  not  too  large,  the  solution  has  the  form 
F(2,t)e‘“‘,  where  F  can  be  chosen  real,  and  where  the  rate  of  variation  of  F  is  slow 
compared  with  u,  that  is  {lfF)dF/dt  u>.  Then  we  have  two  equations  {F  being  real): 


d^F 

dz^ 


—u^dfiF  —  n 


dt 


dF 

e'^-\-uj€"F  =  0 
ut 


(2a) 


(2b) 


where  the  second  time  derivative  of  F  has  been  neglected.  Substituting  (2b)  into  (2a) 
gives  an  equation  for  the  z-variation  of  F 


d^F 

dz^ 


+  w 


F  =  0. 


However,  (2b)  has  a  solution  that  decays  or  grows  indefinitely  depending  on  the  sign  of 
e".  To  determine  a  steady  state,  proceed  as  follows:  Take  the  complete  steady  state  of  the 
field  to  have  the  form 


E  =  Fi  sin  +  Fj  sin 

and  from  now  on  measure  distance  in  units  of  c/x,  defining  C  =  s’r/a.  The  dielectric 
constant  depends  on  E.  Evidently  a  steady  state  would  be  possible  if  d'{E)  •  F  =  0,  for 
each  of  thq  two  modes.  Two  equations  for  the  amplitudes  Fx  and  F2  are  then  obtained  by 
using  trigonometric  identities  to  segregate  out  terms  that  go  like  sin(T2/a)  exp  and  like 
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sin(2xz/c)exp(2iw<)  and  equating  them  to  zero  separately.  The  first  step  is  to  express  e" 
in  terms  of  a  solution  of  Eq.  (1). 

Equation  (1),  in  the  present  one-dimensional  case  is,  in  the  steady  state, 

-J'  =  —  +  iy^-aF^n 

— n  /Pn 

=  -^  +  D'-^-a[F^sm^C  +  F^sm^2C]n  (3) 

where  D'  =  Dv^/a'^  and  J'  =  J/6. 

From  this,  and  from  the  fact  that  sin^^  =  (1/2)(1  —  cos  25),  it  is  apparent  that  a 
consistent  steady  state  is  reached  if  it  is  assumed  that  n  is  a  Fourier  cosine  series  of  the 
variable  2rr/a: 

n  =  no  +  ni  cos  2C  +  n2  cos  4C  +  . . . 

We  substitute  this  into  (3),  use  sum  and  difference  formulae  for  the  cosine  functions,  and 
equate  coefficients  of  the  resulting  trigonometric  functions.  This  gives 

J't  =  [1-f  Ta{F^  -b  F2^)]no  -  ^raFlnx  —  ^rQF’2^n2 
0  =  —TaF^rio  -f  ^1  -4- 4tD'  -f  ra  +  ^-^2^ j  ~ 

0  =  -raF^no  —  ^raFf  Ui  +  [l  +  16ri>'  -f  ra  ^^Fi*  +  F|^  j  712-  (4) 

In  this  process,  the  sum  formulae  also  generate  cos  6C  terms,  but  these  are  discarded. 
They  could  be  taken  into  account  by  including  them  in  the  ansatz  for  n.  Continuing  in 
this  manner,  one  would  eventually  establish  the  full  Fourier  series  for  n.  Hopefiilly  this 
converges  sufficiently  rapidly  so  that  the  first  few  terms  suffice. 

Equations  (4)  are  three  simultaneous  linear  equations  for  no, n2,  which  may  be  solved 
in  terms  of  Fi  and  F?.  If  the  nonradiative  recombination  rate  is  fast  compared  with  the 
induced  rate,  then  raF^  is  small,  amd  an  approximate  solution  to  first  order  in  raF^  should 
suffice: 

no  =  JV(1  —  A  —  B) 

_  Ano 

“  {1+4tD') 
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”2 


Bno 

(1  +  16ri>') 


(5) 


ff2  —  [■A  + 

where 

_  _  [1-70 

are  assumed  positive  (i.e.  gain  rather  than  loss)  and 

9i  =  l  +  i(l  +  4ri?0'\  ?2  =  l  +  i(l  +  16rr>')“'  • 

Equation  (6)  may  be  solved  for  A  and  B: 

4  =  gzffi  —  92 
9i92  -  1 

Q  _  gig2  —  ffi 
gift  - 1  ’ 


(6) 


o 


This  solution,  based  on  the  approximate  solution  (5)  of  Eq.  (4),  breaks  down  when  the  two 
inequalities,  qigi  >  g2  and  q2g2  >  gl.  are  inconsistent.  Evidently,  they  are  consistent  as  long  as 
qiq2  >  1.  For  example,  with  A%D" «  1,  we  have  both  q’s  approximately  equal  to  two,  so  that  the 
inequalities  are  satisfied. 

The  energy  density  at  position  C  is 

W  —  A  sin^  C  +  -^  sin^  2^. 

Figure  2(a)  shows  W  for  equal  gains  in  the  two  modes,  as  a  function  of  C  and  of  r  =  4rZ)', 
which  is  effectively  the  ratio  of  recombination  time  to  the  time  needed  to  diffuse  through 
about  one  wavelength  of  the  lower  mode.  Figures  3  cind  4  show  the  results  for  unequad 
gains,  as  indicated  in  the  captions.  Figures  5  through  7  show  corresponding  carrier  density 
variations. 


III.  Generalization  to  More  Realistic  Structures 


Consider  the  geometry  of  Fig.  8  representing  a  one-stripe  laser  diode.  A  very  similair 
procedure  may  be  used  in  this  case  in  which  the  current  J,  and  therefore  the  dielectric 
constant  is  a  function  of  i.  We  change  that  variable  to  f  =  xr/a.  The  diffusion  equation 
is  then 

-Tj'(f)  =  -n  +  iJ'(^+0)-2aT^F„  (S) 

and  we  write 

n  =  no(0  +  ni(^)sinC  +  n2(Osin2C  +  ....  (9) 


The  equation  for  the  E-field  is 

^  d^E  d^E  f  d^E 

Here  Cq  is  the  background  dielectric  constant,  and 

C)  =  (^)  “  *Ti.2)n(^,C)  +  t‘7o] 


(10) 
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Fig.  2  Total  energy  density  in  the  two  lowest  modes  of  a  parallel  plate  cavity,  a  fimction 
of  distance  z  across  the  cavity  and  of  r,  the  ratio  of  recombination  time  to  the  time 
needed  to  diffuse  through  roughly  one  wavelength.  Here  the  two  modes  have  equal 
gain,  and  are  always  in  the  shaded  region  of  Fig.  1  no  matter  how  large  r:  (a)  view 
from  the  side  r  <  0.  (b)  view  from  large  r. 
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Fig.  7  Same  as  Fig.  10,  but  with  51  =  1,  5'2  =  1-2.  (a)  and  (b)  are  front  and  rear  views. 
Note  that  the  majclmum  r  for  which  n  is  still  positive  is  larger  (up  to  r  =  1.5)  than 
the  maximum  (.375)  allowed  to  ensure  positive  squared  field  amplitudes.  The  latter 
figure  of  course  governs. 
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I  for  a  one-stripe  laser  diode. 
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is  the  locaHzed  deviation  from  Co  induced  by  the  carrier  density.  For  the  steady  state,  we 
assume  a  solution  of  (10)  in  the  form 

E  =  €^^71(0  sin  C  +  e'^72(0  sin2C, 

where  wi  and  wiU  differ  sUghtly  .from  wq  =  CKHa^fe^  and  Tmo-  In  the  Hnear  limit, 
i.e  well  below  threshold,  the  functions  h{i)  and  S^{i)  describe  the  transverse  variations 
of  two  characteristic  modes  of  the  structure  that  go  with  the  z  variations  sin^  and  sin2e, 
respectively.  Here  we  shall  discuss  only  two  possibilities:  (a)  fx  =  lowest  mode  even  m 
/2  =  lowest  mode  odd  in  (b)  both  fx  and  h  lowest  even  modes.  In  exactly  the  same  way 

as  for  the  parallel  plate  cavity,  we  find 

nE  =  (no  -  ini)  sin  C  +  (^o  -  1^2)  sin  2C  +  -  •  • 

where  f\  and  /2  satisfy 


Also,  substituting  (9)  in  (8),  we  find 

rJ'd)  =  l+aT(/f  +  /|)-rD'^jn„-iar/;n,-iaT/|r.j 
0  =  — ar/i  no  +  ^1  +  (/i  +  9/2)  + 

0  =  -arflno  -  iar/f  ni  +  +  orr  (-/i  +  /|)  +  -  tD  nj.  (12) 

As  before,  we  discuss  these  equations  for  nonradiative  recombination  time  short  compared 
with  radiative  recombination  time,  i.e.  ar/?^  «  1.  Then  the  field  terms  in  Eq.  (12)  can 
be  dropped,  except  where  they  compete  with  the  diffusion  term,  i.e.  inside  the  square 
parentheses.  There  they  could  be  important  for  multistripe  lasers.  If  we  were  to  drop  all 

field  dependent  terms,  the  equation  for  no  would  be 

_,,5^no 

T  J  {0  =  no  —  tD 


whose  solution  is 


and  the  exponential  tends  to  fuzz  the  edges  of  the  current  profile.  So  if  two  stripes  are  less 
than  a  distance  VWrx  a  wavelength  apart,  they  would  tend  to  act  like  a  single  stripe.  If 
the  field  terms  are  retained  in  the  parenthesis,  they  will  have  (very  roughly)  the  effect  of 
replacing  r  by  r/[l  +  ra  x  an  average  of  (/f  +  /|)].  It  is  plausible  to  suppose  that  this  is 
the  most  sensitive  spot  at  which  the  fields  affect  the  carrier  concentration,  and  so  we  shall 
disregard  terms  in  higher  powers  of  rap  everywhere  except  this  place. 

We  begin  by  considering  this  problem  for  D'  =  0.  Then  in  the  same  approximation  as 
for  the  parallel  plate  problem 

no  =  rJ'iOil  -  c^rip  +  p)] 
ni  =  aTpno  =  aT^J'{^)p 

n.  =  ar/|no  =  QrrV'(0/|  (13) 

to  first  order  in  arp.  Equations  (11)  now  become 

_  £7i)tJ'(0  (l  -  |ar/f  -  /i  = 

+  [(^  -  i7i)Tj'(0  (l  -  ar/2  -  lar/l)  A  =  <^2/2  (14) 


where  ctu  =  —(1  +  f7o)(‘^i.2/wo)^. 


We  “solve”  these  in  a  manner  similar  to  that  for  the 


parallel  plates,  except  that  we  have  here  no  comparably  convenient  sum  and  difference  formulae 
for  the  eigenfunctions  of  the  7'(^)  potential.  Let  ^0.1(2;)  be  the  lowest  and  next  to  lowest 
eigenfunctions  of  that  potential,  satisfying 

^nj'(^)uui  =  Ao.iV>o,i 

and  going  to  zero  at  ^  =  ioo. 

Suppose  that  J'  is  an  even  function  of  Then  the  eigenfunctions  are  either  even  or 
odd  about  ^  =  0,  and  the  lowest  one  is  even,  the  next  to  lowest  is  odd.  Assume,  with  case 
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/i(0  =  P^oiO  +  terms  in  higher  even  eigenfunctions 
f2i0  —  +  terms  in  higher  odd  eigenfunctions 


(15) 


where  p  and  q  are  constants.  Next  we  expand  /f  etc.  in  a  series  of  even  eigenfunctions, 
and  keep  only  the  first  term.  K  the  eigenfunctions  are  normadized 

noiMo?  =  +■■■■, 

AiHMO?  =  +  ■■■■,  A0{M()Vf2(0  =  p‘ih2M0,  (16) 


where 


^11  “  •4n  J ^21  —  Jm  J  —  ^12>  ^22  —  Jm  j 

and  where  J'{0  =  is  the  maximum  in  the  current  distribution.  We  further 


write 


‘^'(f)/i,2(^)  =  P^i.20i(O.  where  ^1,2  =  J 

The  condition  for  a  steady  state  is  that  the  imaginary  part  of  the  potential  term  should 
vanish  in  both  Eqs.  (11).  This  gives 


+  ^12^1  =  9i  ~ 
^Ai2-P  +  ^^22^1  =  92  — 


7iaTj; 


r  = 


= 


where  P  =  1/2  arp^  etc.  These  may  be  solved  for  P  and  Q  and  presumably  give  very 
simileir  behavior  with  respect  to  the  transverse  dimension  that  we  had  for  the  longitudinal 
dimension  in  the  parallel  plate  case.  We  have 


P  = 


Q  = 


(^2^22^1  —  A12G2) 
A 

(1^11^2  -  Ai2<?i) 


with  A  =  (9/4)^11^22  —  A?2  which  is  always  positive. 
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We  now  include  diffusion.  We  note  that  the  operator  does  not  change  the 

parity  of  a  wavefunction.  We  may  consistently  assume  that  no,  ni,  expand  in  even 
eigenfunctions,  and  we  assume  that  the  lowest  of  these,  ^o,  dominates.  (This  is  a  little 
different  from  the  C-dependence  :  In  the  C  direction  the  densities  persist  right  to  the  end 
plates,  whereas  in  directions  the  n’s  fall  off  to  zero,  like  J\  So  in  the  C  direction  we 
needed  functions  that  were  fimte  at  the  endpoints^  the  cosines  in  this  case.)  This  means 
that  a  term  such  as  D'<Pno,i/d^^  may  be  replaced  by  uq^x^D'  J  =  -25no.i,2, 

where  D  is  positive  because  the  curvature  of  ipx  is  negative  The  only  effect  this  would  have 
on  Eqs.  (13)  is  to  replace  r  by  r/(l  +r!D’)  everywhere.  As  far  as  diffusion  in  the  ^direction 
is  concerned,  the  solution  (IS)  therefore  remains  intact,  except  for  that  replacement  in  the 
definitions  of  P,  Q,  and  the  G’s.  However,  diffusion  in  the  C  direction  is  another  matter. 
Just  as  in  the  case  of  the  parallel  plates,  it  has  the  effect  of  replacing  the  factor  3/2  in  the 
first  of  Eq.  (18)  by  52,  and  in  the  second  of  (18)  by  qi. 

To  plot  the  energy  density  in  the  cavity  as  a  function  of  ^  and  ^  we  need  the  two  lowest 
eigenfunctions.  Suppose  that  J'{^)  has  a  sech^^  variation.  This  has  the  advantage  that  the 

exact  eigenfunctions  are  known.7  We  write  iO)f/co^)zf  (x)  =  r7i,2sech2|.  The  lowest  eigen¬ 
value  is 


and  the  corresponding  (unnormalized)  eigenfunction  is 


The  next  one  has 


with 


V’o  =  (sechO^®- 


3 

2 


rpi  =  sinh^(sechf)*’'^^ 

=  sinh^i/;o 

it  being  assumed  that  q  is  large  enough  to  support  this  mode.  For  convenience,  chose 
V  —  15/4.  Then  Aq  =  3/2  and  Aj  =  1/2.  Then  numerical  evaluation  gives  hn  =  5t/Z2, 
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hi2  ~  ^22  =  ’r/32,  gi  =  2/3,  52  =  1/3-  The  local  energy  density  is 

W  =  Psech^(f  sin^  C  +  Q  sinh^  ^sech^^  sin^  2^.  (19) 

Figures  9  through  11  show  various  cases,  with  fi  and  /2  the  lowest  modes  of  opposite  parity. 
Figures  12  and  13  refer  to  the  ceise  of  fi  and  /2  both  of  even  parity. 

When  r  becomes  too  large  (which  effectively  causes  the  field  induced  recombination 
rate  to  dominate),  this  simple  theory  dependent  on  the  convergence  of  a  mode  expansion 
breaks  down.  This  does  not  necessarily  mean  that  no  steady  state  exists.  If  it  does  exist, 
nonperturbative  solutions  for  the  steady  state  must  be  sought. 
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Fig.  9  Energy  distribution  with  /i  even,  /j  odd,  in  a  single  stripe  structure.  Energy 
density  versus  x  (dimension  transverse  to  the  stripe),  and  z.  Here  Gi  =  SGj  and 
r  =  0.  (In  the  text,  this  implies  that  the  7’s  must  be  different  for  the  two  modes, 
otherwise  the  relation  between  Gi  and  G2  would  be  predetermined.  Thus  for  70  =  0 
in  both  modes  we  would  have  Gi  =  2G2). 


Fig.  10  As  in  Fig.  5,  but  with  Gi  =  2G2  and  r  =  0. 
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Fig.  11  As  in  Fig.  5,  but  with  G\  =  G2>  and  r  =  5.  The  reason  why  the  holebuming 
here  is  generally  more  marked  them  in  the  parallel  plate  cavity  is  that  the  important 
overlap  integrals  of  the  two  modes  with  each  other  and  with  the  current  distribution 
are  different. 


0.1 
0.0 

0. 

0. 

-4 


Fig.  12  /i  and  h  both  the  lowest  even  modes  for  the  single  stripe  case,  for  equal  gains 
and  r  =  0. 
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Fig.  13  Same  as  Fig.  12,  but  with  r  —  5. 
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TECHNOLOGY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer"  for  national  security 
programs,  specializing  in  advanced  military  space  systems.  The  Corporation’s  Technology 
Operations  supports  the  effective  and  timely  development  and  operation  of  national  security 
systems  through  scientific  research  and  the  application  of  advanced  technology.  Vital  to  the 
success  of  the  Corporation  is  the  technical  staffs  wide-ranging  expertise  and  its  ability  to  stay 
abreast  of  new  technological  developments  and  program  support  issues  associated  with  rapidly 
evolving  space  systems.  Contributing  capabilities  are  provided  by  these  individual  Technology 
Centers: 

Electronics  Technology  Center:  Microelectronics,  solid-state  device  physics, 

VLSI  reliability,  compound  semiconductors,  radiation  hardening,  data  storage 
technologies,  infrared  detector  devices  and  testing;  electro-optics,  quantum  electronics, 
solid-state  lasers,  optical  propagation  and  communications;  cw  and  pulsed  chemical 
laser  development,  optical  resonators,  beam  control,  atmospheric  propagation,  and 
laser  effects  and  countermeasures;  atomic  frequency  standards,  applied  laser 
spectroscopy,  laser  chemistry,  laser  optoelectronics,  phase  conjugation  and  coherent 
imaging,  solar  cell  physics,  battery  electrochemistry,  battery  testing  and  evaluation. 

Mechanics  and  Materials  Technology  Center:  Evaluation  and 
characterization  of  new  materials:  metals,  alloys,  ceramics,  polymers  and  their 
composites,  and  new  forms  of  carbon;  development  and  analysis  of  thin  films  and 
deposition  techniques;  nondestructive  evaluation,  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  development  and  evaluation  of 
hardened  components;  analysis  and  evaluation  of  materials  at  cryogenic  and  elevated 
temperatures;  launch  vehicle  and  reentry  fluid  mechanics,  heat  transfer  and  flight 
dynamics;  chemical  and  electric  propulsion;  spacecraft  structural  mechanics,  spacecraft 
survivability  and  vulnerability  assessment;  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  lubrication  and 
surface  phenomena. 

Space  and  Environment  Technology  Center:  Magnetospheric,  auroral  and 
cosmic  ray  physics,  wave-particle  interactions,  magnetospheric  plasma  waves; 
atmospheric  and  ionospheric  physics,  density  and  composition  of  the  upper 
atmosphere,  remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared 
astronomy,  infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere;  effects 
of  electromagnetic  and  particulate  radiations  on  space  systems;  space  instrumentation; 
propellant  chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
atmospheric  chemical  reactions,  atmospheric  optics,  light  scattering,  state- specific 
chemical  reactions  and  radiative  signatures  of  missile  plumes,  and  sensor  out-of-field- 
of-view  rejection. 


